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ABSTRACT 


Whistler  waves  are  a  type  of  very  low  frequency  (VLF)  radiation  which  propagate 
through  the  earth’s  magnetosphere.  This  phenomenon  can  occur  as  a  result  of  a 
lightning  discharge  and  can  be  produced  using  ground-based  VLF  transmitters.  It  is 
theorized  that  these  signals  travel  through  ducts  centered  no  geomagnetic  field  lines. 
While  passing  through  these  ducts  the  waves  are  amplified  through  resonant 
interaction  with  electrons  in  the  radiation  belts  which  are  following  helical  paths  around 
magnetic  field  lines.  A  description  of  whistlers  and  the  related  topic  of  VLF  emissions 
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is  presented  in  Chapter  II,  along  with  a  description  of  how  these  effects  change  the 
properties  of  the  ionosphere. 

Starting  with  the  Lorentz  force  equation,  the  equations  of  motion  for  electrons  in 
the  magnetosphere  are  developed  in  Chapter  III.  These  equations  are  numerically 
integrated  along  with  a  wave  equation  developed  in  Chapter  IV.  The  results  of  the 
simulation  are  presented  in  Chapter  V  for  the  idealized  case  of  a  monoenergetic  beam 
of  electrons.  High  and  low  gains  in  the  wave  amplitude  are  observed  for  both  strong 
and  weak  initial  VLF  fields.  The  simulation  is  also  run  using  distributions  to  represent 
the  initial  energies  of  the  electrons  in  the  radiation  belts.  Chapter  VI  presents  some 
possible  ways  to  make  the  simulation  more  realistic  along  with  a  summary  of  the 
similarities  between  the  theory  presented  and  Free  Electron  Laser  theory. 
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I.  INTRODUCTION 


The  study  of  "whistlers"  dates  back  to  the  earliest  days  of  radio 
communications.  The  first  reported  survey  of  these  distinct  musical  type  signals 
came  about  as  a  result  of  German  attempts  to  listen  in  on  Allied  telephone 
conversations  during  World  War  I  [1].  From  these  beginnings,  has  developed  a 
theory  of  how  whistlers  are  produced  and  how  the  waves  propagate  between 
magnetic  conjugant  points  in  the  Northern  and  Southern  hemisphere.  This 
information  is  reviewed  in  Chapter  II  of  this  thesis.  The  study  of  the  related 
phenomenon  of  VLF  emissions  dates  back  to  the  early  1 950’s  and  a  review  of 
this  theory  is  also  presented  in  Chapter  II.  The  chapter  concludes  with  a  look 
at  how  these  signals  can  limit  the  range  of  VLF  communication  systems  and 
degrade  the  accurracy  of  the  navigation  systems  used  by  the  U.  S.  Navy 

In  Chapter  III  a  description  of  the  magnetosphere  is  provided.  Most  of  the 
early  information  on  this  region  of  space  comes  from  the  study  of  whistler 
spectrograms.  Chapter  III  concludes  with  a  discussion  about  the  electrons  in 
the  magnetosphere  along  with  the  equations  of  motion  used  to  describe  how 
the  electrons  move  in  the  geomagnetic  field. 

A  theory  proposed  to  explain  the  amplification  of  VLF  signals  in  the 
magnetosphere  is  that  of  wave-particle  interactions  between  incoming  VLF 
waves  and  a  stream  of  counter-propagating  electrons.  This  resonant  interaction 
is  thought  to  occur  in  a  region  of  space  within  a  few  degrees  above  and  below 
the  geomagnetic  equator,  and  at  geocentric  distances  of  three  to  six  earth  radii. 
The  rest  of  the  thesis  explores  this  nonlinear  interaction  through  the  use  of 
computer  simulations.  In  Chapter  IV,  background  plasma  in  the  interaction 
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region  is  described  which  leads  to  a  wave  equation  used  to  follow  the 
development  of  the  amplitude  and  phase  of  the  wave  as  it  passes  through  the 
interaction  region.  The  motion  of  the  resonant  electrons  in  the  presence  of  the 
wave  is  described  using  a  set  of  first  order  differential  equations  which  are 
developed  from  the  Lorentz  force  equation. 

The  computer  simulation  numerically  integrates  the  electron  equations  of 
motion  and  the  wave  equation  given  a  set  of  initial  conditions.  In  Chapter  V, 
the  results  of  this  simulation  are  presented  for  a  variety  of  initial  conditions. 
The  first  set  of  simulations  explore  what  happens  for  the  case  of  a 
monoenergetic  beam  of  electrons.  Conditions  needed  to  produce  high  and  low 
gain  in  the  wave  are  explored  as  well  as  strong  and  weak  fields.  The  second 
set  of  runs  uses  initial  conditions  which  approximate  those  found  in  the 
ionospheric  interaction  region.  The  initial  phase  velocities  and  transverse 
velocities  for  the  resonant  electrons  are  determined  using  distribution  functions. 

In  Chapter  VI,  suggestions  for  future  work  to  improve  the  simulation  is 
discussed.  It  turns  out  that  the  theory  presented  in  this  thesis  regarding  wave- 
particle  intercations  has  many  similarities  to  the  theory  used  to  describe  the 
Free  Electron  Laser.  These  similarities  are  discussed  at  the  end  of  Chapter  VI. 
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II.  VLF  RADIATION  IN  THE  MAGNETOSPHERE 


A.  WHISTLERS 

"Whistler  waves"  are  one  type  of  very  low  frequency  (VLF)  electromagnetic 
radiation  which  propagate  through  the  earth's  magnetosphere  and  are  detected 
at  magnetically  conjugate  points  on  the  earth’s  surface.  The  name  "whistler" 
provides  the  best  description  of  what  these  emissions  sound  like;  a  high  pitched 
tone  (in  the  kilohertz  range)  which  sweeps  down  in  frequency  over  the  period  of 
about  one  second.  Whistlers  can  be  heard  by  the  human  ear  using  very  simple 
equipment.  As  a  result,  reports  of  their  occurrence  date  back  to  the  earliest 
radio  experiments.  Two  metal  probes  placed  in  the  ground  50  feet  apart  and 
connected  to  a  high  gain  amplifier  makes  an  excellent  whistler  receiver.  They 
can  also  be  heard  over  telephones  connected  to  long  telephone  lines.  The  rate 
of  whistler  activity  varies  with  the  time  of  day,  magnetic  activity  and  location  on 
the  earth.  They  are  more  common  at  night  and  frequently  occur  in  conjunction 
with  lightning  storms.  Mid-latitude  regions  report  the  highest  rate  of  whistler 
activity,  while  lower  rates  are  found  in  the  higher  latitudes.  Whistlers  have  not 
been  observed  at  the  equator  or  at  the  poles.  The  latitudes  where  the  highest 
rates  are  reported  seem  to  move  toward  the  equator  during  periods  of 
increased  magnetic  activity.[1] 

The  close  association  between  lightning  discharges  and  the  detection  of 
whistlers  has  lead  to  the  following  explanation  for  the  occurrence  of  these 
waves.  A  lightning  discharge  releases  radio  waves  with  a  wide-band  of 
frequencies.  Some  of  these  waves  reflect  off  the  ionosphere,  at  altitudes 
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between  90  and  200  kilometers  (km)  and  can  travel  over  the  horizon  to  ranges 
of  thousands  of  kilometers.  It  is  postulated  that  some  of  the  radiation  from  the 
lightning  penetrates  the  ionosphere  and  enters  ducts  which  are  centered  on 
geomagnetic  field  lines.  The  waves  follow  these  field-aligned  ducts  from  one 
hemisphere  to  the  magnetic  conjugate  point  in  the  other  hemisphere.  The 
waves  pass  through  the  ionosphere  and  are  detected  as  whistlers.  The 
distinctive  falling  pitch  indicates  that  the  path  taken  by  the  whistlers  is  highly 
dispersive.  Higher  frequency  components  of  the  wave  travel  faster  through  the 
ducts  than  those  with  a  lower  frequency,  arriving  at  the  receiver  first.  A  portion 
of  the  wave  can  also  reflect  off  the  top  of  the  ionosphere,  returning  to  the 
originating  hemisphere  along  the  same  duct.  This  "two-hop"  whistler  is  then 
detected  at  a  ground  receiver.  This  reflection  between  magnetic  conjugate 
points  can  occur  many  times  and  produce  a  number  of  echos.[2]  In  addition  to 
lightning,  VLF  radiation  from  a  nuclear  explosion  and  man-made  Morse  code 
dashes  from  VLF  transmitters  can  all  travel  in  this  same  whistler  mode  [1]. 

The  ducts  which  whistlers  travel  through  are  believed  to  work  on  the  same 
principles  as  a  fiber  optic  cable.  There  is  a  region  around  the  magnetic  field 
line  of  increased  electron  density.  This  results  in  a  higher  index  of  refraction 
and  the  electromagnetic  waves  tend  to  stay  within  this  region  due  to  total 
internal  reflection,  as  long  as  the  angle  of  incidence  between  the  wave  and  the 
wall  of  the  duct  is  small.  Extensive  study  of  whistler  data  have  shown  that  the 
time  between  successive  whistler  echos  at  one  receiver  is  always  an  integral 
number  of  seconds.  In  addition,  the  dispersion  properties  of  two  whistlers  do 
not  change  over  the  course  of  many  minutes  indicating  that  the  whistlers  travel 
the  same  path.  These  results  reinforce  the  theory  that  the  waves  follow  ducts 
which  are  aligned  with  magnetic  field  lines. [3] 
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A  convenient  way  to  display  whistlers  is  with  a  frequency  verses  time 
spectrogram  shown  in  Figure  2-1.  The  small  arrow  indicates  reception  of  the 
impulse  which  causes  the  whistler,  possibly  a  stroke  of  lightning  in  the  case  of 
naturally  occurring  whistlers.  The  dark,  curved  trace  indicates  reception  of  the 


Figure  2-1 :  Whistler  spectrogram  on  a  frequency  versus  time  scale  [1]. 

whistler  signal.  The  darker  the  line  the  higher  the  intensity  received.  The  time 
delay  between  the  impulse  and  reception  of  the  whistler  is  proportional  to  the 
length  of  the  whistler’s  path  of  propagation.  In  this  case  the  first  trace  is  a 
one-hop  whistler.  It  has  made  one  pass  through  the  duct.  The  smooth 
decrease  in  frequency  is  clearly  evident.  Also  shown  in  Figure  2-1  is  an  echo 
or  "three-hop"  whistler  which  has  made  three  trips  through  the  same  field 
aligned  duct  before  passing  through  the  ionosphere  and  being  picked  up  by  the 
receiver.  The  rate  of  frequency  decrease  in  this  and  succeeding  whistler  echos 
is  not  as  great  as  in  the  initial  whistler  trace.  Each  successive  echo  makes  one 
more  pass  through  the  dispersive  medium  in  the  duct  then  the  previous  echo, 
stretching  out  the  spectrogram  for  each  successive  echo.[1] 
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The  initiating  lightning  impulse  does  not  sweep  in  frequency  because  it 
travels  by  a  direct  path  in  the  earth  ionosphere  waveguide  to  the  receiver.  Very 
little  dispertion  takes  place  below  the  ionosphere  and  the  travel  time  for  the 
impulse  is  a  fraction  of  a  second  resulting  in  a  spike  made  up  of  many 
frequencies  and  very  short  in  duration.[1] 

Figure  2-2  shows  a  group  of  whistlers  packed  very  close  together.  These 
are  an  examples  of  "multipath  whistlers".  Some  of  the  energy  coming  from  a 
particular  lightning  stroke  can  enter  adjacent  ducts,  each  with  a  slightly  different 
path  length.  As  a  result  the  waves  have  different  travel  times  and  when 
received  appear  to  overlap  as  shown  in  Figure  2-2. [4] 


Figure  2-2:  Example  of  multipath  whistlers  [1]. 

t{( a)  =  f-^-  2-1 

N. 

where  vg  is  the  whistler  group  velocity  and  ds  is  an  increment  of  path  length. 
This  travel  time  is  a  function  of  frequency  (©)  because  plasma  in  the  duct 
provides  a  dispersive  medium.  Assuming  that  the  group  velocity  is  constant 
over  the  total  path,  it  can  be  taken  out  of  the  integral  and  written  as 
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v„  = 


2-2 


where  c  is  the  speed  of  light  in  a  vacuum  and  ng  is  the  group  index  of  refraction 
which  is  defined  as  ng  =  n  +  (o  dn/da>.  The  index  of  refraction,  n,  is  also  a 
function  of  frequency  and  is  known  to  have  the  form  [3] 


n  s  o-3 

V(o(ft-w)  ’ 

where  the  plasma  frequency,  ©c,  and  the  electron  cyclotron  frequency,  ft,  are 
defined  as 


(Or  = 


4ite2p 


m 


and  ft  = 


eB 

me 
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where  e  is  the  electron  charge  magnitude,  B  is  the  magnitude  of  the  magnetic 
field,  m  is  the  electron  mass  and  p  is  the  electron  density.  The  approximation 
stated  in  equation  2-3  is  accurate  for  the  case  where  <ac  »  to  and  <oc  »  ft. 
This  situation  exists  for  the  region  where  whistlers  propagate.  Assuming  that  the 
waves  follow  a  particular  geomagnetic  line  of  force  and  that  the  earth’s 
magnetic  field  is  a  dipole,  the  path  length  can  be  determined  [1].  The  resulting 
equation  for  travel  time  is 


f(a>)  = 


/  <ocn 

2c  o)1/2 


(ft-co)_3/z  . 
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where  /  is  the  total  path  length.[5] 

Figure  2-3  is  a  plot  of  f(<o).  In  this  graph,  frequency  is  normalized  to  the 
electron  cyclotron  frequency  and  time  is  given  in  seconds.  Also  indicated  is  the 
frequency  of  minimum  time  delay,  know  as  the  "nose  frequency"  It  can  be 
shown  using  the  equation  for  f(a>)  that  a  relationship  exists  between  <*>„  and  ft: 


7 


<d„  =  0.25ft  . 
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Substituting  this  value  back  into  the  equation  for  f(<o)  results  in  an  expression 
for  the  nose  time  delay  tn : 


8V3  /  (0C 
tn  =  3  c  ft 
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This  derivation  results  in  a  way  to  find  ft  using  the  value  for  (o„  which  is  read 
from  the  nose  whistler  spectrogram.  Then  for  a  given  path  length  /,  the 
plasma  frequency  ac  can  be  determined  using  equation  2-7.  This  analysis 
assumes  the  plasma  in  the  magnetosphere  is  uniformly  distributed.[5] 


Figure  2-3:  Plot  of  t( to)  [5]. 

The  real  magnetosphere  does  not  completely  comply  with  the  assumptions 
made  in  the  preceding  analysis.  Extensive  study  of  nose  whistlers,  such  as 
pictured  in  Figure  2-4,  has  lead  to  an  accurate  relationship  between  ©„  and  ft: 
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Figure  2-4:  Nose  whistler  spectrogram  [1]. 


B.  VLF  EMISSIONS 

A  second  type  of  wave  that  propagates  through  the  magnetosphere  is 
known  as  a  VLF  emission.  These  waves  fall  into  two  major  categories,  "hiss" 
shown  in  Figure  2-5  and  discrete  emissions  also  shown  in  Figure  2-5  and  which 
go  by  the  names  "risers",  "fallers"  and  "hooks"  [1j.  Hiss  is  broad  band  low 
frequency  radiation  which  has  a  definite  cut-off  frequency.  The  various  types  of 
discrete  emissions  are  narrow  band  and  can  sweep  through  a  range  of 
frequencies  over  the  period  of  a  second.  The  same  equipment  is  used  in  the 
study  of  VLF  emissions  as  is  used  to  detect  whistlers.  VLF  emissions  can 
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occur  in  conjunction  with  strong  whistlers.  A  whistler  can  trigger  a  discrete 
emission  such  as  in  Figure  2-6,  where  the  emission  begins  at  the  low  frequency 
tail  of  the  whistler.  These  types  of  whistler  triggered  emissions  are  more 
frequent  during  periods  of  high  magnetic  activity  and  are  limited  to  the  higher 
latitudes.  In  the  polar  regions,  VLF  emissions  can  significantly  degrade  the 
performance  of  communications  systems.  They  exhibit  the  same  level  of 
intensity  as  whistlers  and  are  also  more  common  in  middle  to  higher  latitudes. 
There  is  a  close  relationship  in  form  and  time  between  emissions  detected  at 
magnetic  conjugate  points  on  the  earth’s  surface.  This  observation  has  led  to 
the  theory  that  the  waves  travel  through  the  magnetosphere  in  the  same 
manner  described  for  whistlers.  This  is  known  as  the  whistler-mode  of 
propagation.  There  is  also  strong  association  between  auroral  activity  and  an 
increase  in  VLF  emission  rates.[1] 

The  sources  of  VLF  emissions  are  not  completely  understood.  They 
probably  originate  in  the  magnetosphere  and  are  the  result  of  the  interaction 
between  waves  propagating  in  the  whistler-mode  and  high  energy  electrons 
(keV)  located  in  the  radiation  belts  which  spiral  around  geomagnetic  field  lines. 
In  addition  to  whistler  triggered  emissions  it  has  been  found  that  ground  based 
VLF  transmitters  can  produce  triggered  emissions  as  shown  in  Figure  2-7. 
Extensive  research  using  these  man-made  triggered  emissions  has  led  to  the 
following  list  of  characteristics: 

1 .  A  minimum  pulse  length  is  needed  to  produce  a  triggered  emission.  Only 

Morse  code  dashes  of  length  100-150  ms  result  in  an  emission,  Morse 

code  dots  of  length  50ms  do  not  trigger  an  emission. 
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Figure  2-5:  Examples  of  VLF  emiss 
"risers"  and  "hooks"  [1]. 
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tellers", 


2.  The  original  VLF  signal  shows  growth  in  intensity  on  the  order  of  25-30  dB 
prior  to  the  triggered  emission.  The  triggered  emission  starts  once  the 
initiating  signal  stops.  Observations  of  triggered  emissions  caused  by 
whistlers  shows  that  the  triggered  signal  begins  toward  the  end  of  the 
whistler  pulse  in  the  range  where  the  rate  of  frequency  decrease  is  a 
minimum. 

fcC  _ 


Figure  2-6:  Example  of  a  whistler  triggered  VLF  emission  in  this  case 
a  hook  [1]. 


3.  The  triggered  signal  begins  at  (or  within  200HZ  of)  the  frequency  of  the 
initiating  pulse  and  sweeps  up  or  down  in  frequency  in  a  fraction  of  a 
second.  The  longer  the  initial  pulse  the  more  likely  that  a  riser  will  result. 
The  emission  in  extremely  narrow  band  (=100Hz)  and  the  sweeping 
frequency  change  can  be  as  much  as  50%.  The  triggered  signal  also 
tends  to  end  abruptly. 

4.  There  are  occasions  when  triggered  emissions  do  not  occur  when  a  Morse 
code  dash  is  transmitted.  The  reason  for  this  is  not  known  but  the  gain 
and  frequency  variation  described  above,  when  a  triggered  emission  does 
occur,  "are  relatively  insensitive"  to  the  level  of  geomagnetic  activity.[7] 
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D.  EFFECTS  OF  VLF  EMISSIONS 


The  interaction  which  causes  VLF  emissions  is  also  considered  to  be  the 
mechanism  that  causes  high  energy  electrons  to  precipitate  out  of  the  earth’s 
radiation  belts.  The  electron  spiral  motion  around  the  magnetic  field  lines 
couples  with  VLF  waves  and  amplifies  them.  This  process  reduces  the  energy 


Figure  2-7:  A  constant  14.7kHz  signal  and  the  resulting  triggered 
emissions.  In  this  case  risers,  falters  and  hooks  are  generated  [1]. 

of  the  electrons  and  causes  them  to  enter  the  ionosphere.  These  high  energy 
electrons  can  increase  the  plasma  frequency  in  the  low  levels  of  the  ionosphere 
directly  through  ionization  or  indirectly  through  the  production  of 
bremmstrahlung  x-rays  or  ultraviolet  emissions.  The  result  of  this  is  the 
formation  of  a  "pimple"  protruding  down  from  the  bottom  boundary  of  the 
ionosphere  as  illustrated  in  Figure  2-8  [9],  VLF  waves  passing  through  the 
earth-ionosphere  wave  guide  can  scatter  off  this  region,  adding  constructively 
and  destructively  to  the  unscattered  wave  and  resulting  in  changes  to  the 
amplitude  and  phase  of  the  original  signal.  This  phenomenon,  known  as  a 
"Trimpi  event"  [10]  is  well  documented  and  can  be  initiated  by  naturally 
occurring  whistlers  which  then  trigger  VLF  emissions  and  by  ground-based  VLF 
transmitters. 
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The  VLF  communication  systems  and  radio  navigation  systems  used  by 
the  U.  S.  Navy  are  susceptible  to  changes  in  the  ionosphere  produced  by  the 
Trimpi  effect.  One  such  system  is  the  Omega  radio  navigation  system  in  which 
a  user  determines  his  distance  from  a  particular  transmitter  by  comparing  the 
phase  difference  between  two  signals  of  the  same  frequency  sent  out  from  two 
different  locations.  There  are  eight  Omega  stations  throughout  the  world  each 


disturbed 

region 


transmitting  the  same  series  of  three  VLF  signals  [11].  Observations  of  signals 
used  in  the  Omega  radio  navigation  system  have  shown  phase  perturbations  on 
the  order  of  a  microsecond.  Assuming  the  signals  propagate  at  approximately 
the  speed  of  light  in  a  vacuum  (3xlo10cm/s)  the  resulting  uncertainty  in 
source-to-receiver  distance  would  be  =  300  meters  for  a  one  microsecond 
phase  difference.  This  amount  does  not  significantly  impact  on  the  stated  fix 
accuracy  of  the  Omega  system  which  is  usually  1  to  4  km,  but  the  amount  of 
error  introduced  is  important  when  using  Differential  Omega.  Differential 
Omega  increases  the  overall  accuracy  of  a  fix  by  assuming  all  signals  received 
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have  undergone  the  same  amount  of  phase  perturbation,  an  assumption  not 
necessarily  true  because  Trimpi  effects  are  localized  and  can  effect  only  one  of 
the  received  signals.  In  addition,  studies  of  Omega  signals  show  that  the 
amount  of  phase  shift  above  60°  latitude  is  of  the  order  of  10  microseconds. 
This  translates  into  distance  errors  of  3000  meters  which  would  double  the 
uncertainty  of  the  fix.[1 2] 

There  is  great  interest  in  learning  more  about  how  the  reflective  properties 
of  the  ionosphere  can  be  changed.  The  precipitation  of  electrons  as  previously 
described  is  one  method  to  change  these  properties  of  the  ionosphere.  Other 
methods  include  the  direct  heating  of  the  ionosphere  through  the  use  of  ground 
based  high  powered  VLF  transmitters  and  the  introduction  of  chemicals  to 
increase  the  ion  concentration.  One  application  of  this  work  is  development  of 
an  over-the-  horizon  radar  which  would  use  a  patch  of  increased  ionization  in 
the  atmosphere  as  a  mirror  [13].  Increased  detection  range  and  reduced 
chance  of  counter-detection  are  two  features  making  such  a  system  extremely 
desirable.  A  better  understanding  of  the  wave-particle  interaction  which  results 
in  electron  precipitation  is  a  step  on  the  road  to  the  development  of  this  type  of 
system. 
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III.  THE  ENVIRONMENT  AND  EQUATIONS  OF  MOTION 


A.  THE  MAGNETOSPHERE 

A  good  approximation  for  the  earth’s  geomagnetic  field  is  that  of  a  dipole 
field  located  within  the  earth.  The  magnetopause  is  a  boundary  which 
separates  the  region  where  the  earth’s  magnetic  field  dominates  from  the  region 
where  the  action  of  the  solar  wind  dominates.  Within  the  magnetopause  the 
dipole  field  approximation  is  valid.  Outside  the  magnetosphere,  the  earth’s  field 
is  distorted  due  to  the  solar  wind.  A  bow  shock  forms  on  the  side  facing  the 
sun  as  the  solar  wind  causes  the  magnetosphere  to  be  compressed.  As  the 
solar  wind  passes  over  and  around  the  magnetosphere,  a  tail  is  created  which 
has  been  observed  to  stretch  out  to  200 Re,  where  Rg  is  the  mean  radius  of  the 
earth  (6370  km).  Figure  3-1  is  a  sketch  of  the  magnetosphere  and  shows  how 
the  solar  wind  distorts  the  shape.[5] 

This  thesis  will  concentrate  on  the  region  inside  the  magnetosphere.  This- 
area  includes  the  plasmasphere,  made  up  of  a  "cold"  background  plasma  and 
the  "hot"  plasma  which  makes  up  the  earth’s  radiation  belts.  The  background 
plasma  is  usually  described  as  electrons  and  protons  ranging  in  energy  from  .1 
to  1  eV.  They  are  in  diffusive  equilibrium  and  have  densities  from  104 
particles/cm3  at  ionosphere  heights  (1000km)  to  300-400  particles/cm3  at  the 
boundary  of  the  plasmasphere.  This  boundary,  called  the  plasmapause,  was 
discovered  through  the  study  of  naturally  occurring  whistlers  and  is 
characterized  by  a  significant  decrease  in  cold  plasma  density.  The  position  of 
the  plasmapause  varies  with  solar  activity:  closer  to  the  earth  during  magnetic 
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Figure  3-1  A  schematic  diagram  of  the  Earth's  magnetosphere  in  the 
noon-midnight  plane  [5]. 

storms  (2 Re)  and  farther  away  (7 Ra)  after  many  days  of  little  activity.  The  "hot" 
plasma  includes  electrons  and  protons  having  energies  ranging  from  10  kev  to 
100  Mev  which  are  trapped  by  the  geomagnetic  field.  They  follow  helical  paths 
around  field  lines  and  bounce  between  magnetic  conjugate  points.  Figure  3-2 
is  a  sketch  of  the  magnetosphere  and  includes  examples  of  the  helical  paths 
these  particles  travel.[1 4] 
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Figure  3-2:  The  magnetosphere  including  energetic  particle  trajectory 
[14]. 


B.  THE  MAGNETIC  FIELD 


The  earth’s  magnetic  dipole  field  is  assumed  to  be  embedded  in  the  center 
of  the  earth  with  the  axis  tilted  approximately  11°  to  the  rotation  axis.  The  field 
lines  emerge  from  the  south  geographic  pole  and  reenter  the  earth  at  the  north 
geographic  pole  as  shown  in  Figure  3-3.  It  is  convenient  to  write  the 
components  of  the  geomagnetic  field  in  terms  of  spherical  coordinates, 
substituting  geomagnetic  latitude  X  for  the  co-latitude  e.  The  field  is  symmetric 
in  the  <)>  direction.  The  equations  describing  the  field  components  and  the 
magnitude  of  the  field  are 


Br  =  -B0 


3 

2sinX 


3-1 
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B,  =  Be 


B  =  Bb 


cosX  . 


Vl+3sin2X  . 


3-2 


3-3 


where  B0  is  the  magnitude  of  the  field  on  the  earth’s  surface  at  the 
geomagnetic  equator  (30,293nT)  and  R0  is  the  mean  radius  of  the  earth  (6370 


Figure  3-3:  Dipole  field  used  to  represent  Earth’s  magnetic  field  [14]. 


The  equation  for  a  dipole  field  line  of  force  is  r  =  R^ cos2X  in  which  f?*,  is 
the  radial  distance  to  the  field  line  when  X  equals  zero.  The  Mcllwain  L 
parameter  is  used  to  express  a  given  field  line  in  terms  of  the  radial  distance  at 
which  the  field  line  crosses  the  geomagnetic  equator, 

L  =  .  3-4 

R0 

This  parameter  is  a  convenient  way  to  label  the  field  aligned  ducts  through 
which  whistler-mode  signals  travel.  [1 4] 
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C.  ELECTRON  MOTION  IN  THE  MAGNETOSPHERE 


The  source  of  the  high  energy  electrons  in  the  radiation  belts  is  not  well 
known,  although  one  theory  is  that  they  are  deposited  by  the  solar  wind  in  the 
polar  regions  and  then  propagate  down  the  magnetic  field  lines  [9].  The  motion 
of  these  electrons  is  governed  by  the  Lorentz  equation, 

rwt  =P  =  -e(\7  x  8)  .  3-5 

This  force  is  always  perpendicular  to  ,  the  velocity  vector  of  the  electron  and 
therefore  the  magnetic  field  8  does  no  work  on  the  electron.  Only  the 
component  of  the  velocity  which  is  perpendicular  to  8  is  affected  by  the  force 
resulting  in  uniform  circular  motion  around  the  magnetic  field  line.  The  radius  of 
the  cyclotron  orbit  is 


rc  = 


mVi 

~eB~ 


3-6 


where  e  is  the  charge  of  the  electron,  m  is  the  mass  of  the  electron,  B  is  the 
magnitude  of  the  magnetic  field  and  Vi  is  the  component  of  the  electron 
velocity  perpendicular  to  8.  Assuming  that  8  points  into  the  paper,  then  the 
electron  will  rotate  in  a  clockwise  direction.  When  the  V  has  a  component 
parallel  to  8,  the  electron  follows  a  helical  path  around  the  field  line.  The  angle 
that  makes  with  the  magnetic  field  line  is  known  as  the  pitch  angle  a 


a  =  tan-1 


3-7 


This  angle  will  change  as  Vi  and  V,,  varies.  When  V„  equals  zero,  a  is  90°.  At 
this  point  the  electron  will  stop  its  forward  motion  then  travel  in  the  opposite 
direction  along  the  field  line.  It  is  as  if  the  electron  reflects  off  a  mirror,  while 
the  direction  of  rotation  remains  the  same.[5] 
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The  magnitude  of  the  earth's  magnetic  field,  for  a  given  field  line,  is  at  a 
minimum  at  the  geomagnetic  equator  and  increases  as  you  move  along  the 
field  line  away  from  the  equator.  An  electron  moving  up  a  field  line  from  the 
equator  will  also  see  an  increasing  B  field,  but  the  electron  executes  many 
cyclotron  gyrations  in  a  distance  where  the  amount  of  change  in  the  field  is 
small.  This  is  known  as  an  adiabatic  change  and  it  leads  to  the  concept  of  the 
adiabatic  invariant  which  states  that  the  square  of  the  transverse  momentum  of 
the  electron  divided  by  the  magnitude  of  the  field  (pf  IB)  is  a  constant.  The 
field  does  not  do  work  on  the  electron  hence  the  speed  of  the  electron  must  be 
constant  and  it  can  be  written  in  terms  of  the  velocity  perpendicular  and  parallel 
to  the  geomagnetic  field 

V02  =  V2  +  V*  ,  3-8 

where  V02  =  Vl02  +  V]f)2  is  the  speed  of  the  electron  at  the  geomagnetic  equator. 
Assuming  a  coordinate  system  in  which  the  z  axis  is  aligned  with  the 
geomagnetic  field,  then  the  speed  of  the  electron  at  z= 0  is  V0  and  the 
magnitude  of  the  field  is  B0.  As  a  result  of  the  adiabatic  invariant, 


B(z) 


3-9 


where  B  is  a  function  of  z.  This  relationship  can  then  be  used  to  find  the 
parallel  velocity  of  the  electron, 


3-10 


As  the  field  increases  with  increasing  z  eventually  the  right  side  of  equation  3- 
10  will  equal  zero.  At  this  point  the  electron  reverses  direction  and  heads  back 
toward  the  equator  while  continuing  to  rotate  around  the  field  line  in  the  same 
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direction.  The  equation  of  motion  describing  this  deceleration  can  be  expressed 
to  first  order  in  terms  of  the  divergence  of  the  field  in  the  z  direction  [15] 


iV  -  Vi02  dB(z) 
-V"  2B0  dz 


3-11 


When  the  electron  reverses  direction  the  pitch  angle  a  equals  90°.  Using 
the  adiabatic  invariant  leads  to  the  definition  of  a  mirror  pitch  angle  am : 

sin2am  =  ,  3-12 

where  Bm  is  the  magnitude  of  the  field  when  the  electron  mirrors  [5].  Any 
electron  with  a  pitch  angle  less  than  am  as  it  crosses  the  geomagnetic  equator 
will  not  mirror.  The  solid  angle  subtended  by  am  defines  a  "loss  cone". 
Electrons  with  pitch  angles  within  the  escape  cone  will  mirror  at  a  latitude  where 
the  field  strength  is  greater  than  Bm.  A  higher  magnitude  means  the  mirror 
point  is  closer  to  the  pole  hence,  these  electrons  will  travel  to  a  higher  latitude 
before  they  reflect. 

As  the  mirror  latitude  for  an  electron  increases  the  altitude  above  the  earth 
where  reflection  occurs  decreases.  As  a  result,  electrons  enter  regions  of  the 
upper  ionosphere  of  higher  particle  density  and  the  possibility  of  energy  loss 
through  collision  increases.  If  collisions  do  occur  the  electrons  will  precipitate 
out  of  their  field  aligned  ducts.  The  wave-particle  interaction  which  is  thought  to 
produce  VLF  emissions  takes  energy  away  from  the  electrons  reducing  their 
pitch  angle.  As  a  result,  they  precipitate  out  of  the  radiation  belts  in  much  high 
flux  than  is  usually  observed.  This  phenomenon  is  the  Trimpi  effect  discussed 
in  the  last  chapter.  The  derivation  and  simulations  presented  in  the  following 
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chapters  will  concentrate  on  the  interaction  which  produces  the  increase  in 
electron  precipitation. 

The  equation  of  motion  for  electrons  in  a  dipole  field  can  be  solved 
numerically  using  a  desk  top  computer.  Figure  3-4  shows  the  path  of  a  trapped 
electron  around  the  earth.  In  addition  to  the  trapping  of  electrons,  the 
inhomogeneity  of  the  earth’s  magnetic  field  causes  the  electrons  to  drift  in 
longitude  as  is  clearly  evident  in  the  figure.  The  effect  of  this  drift  on  the  wave- 
particle  interaction  is  not  significant,  and  will  not  be  considered  in  the 
simulations  presented  later. 


Figure  3-4:  Plot  of  path  taken  by  one  electron  in  a  dipole  field. 
Bounce  and  drift  motion  are  present. 
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D.  ELECTRON  EQUATIONS  OF  MOTION  IN  A  MAGNETIC  FIELD 

The  rest  of  this  thesis  will  look  at  the  wave  particle  interaction  which 
amplify  the  incoming  VLF  wave  and  produce  VLF  emissions.  In  this  model,  the 
electrons  which  are  executing  cyclotron  motion  around  the  field  lines  gives  up  a 
small  amount  of  energy  to  the  propagating  VLF  wave  [1].  When  the  Doppler- 
shifted  wave  frequency  seen  by  the  electrons  equals  the  electron  cyclotron 
frequency,  the  resonance  condition  can  be  expressed  as 

k  V,  =  co  -  ft  ,  3-13 

where  k  is  the  wave  number  of  the  VLF  wave.  For  the  case  of  waves 
propagating  in  the  whistler  mode,  the  electron  cyclotron  frequency,  Q,  at  the 
magnetic  equator  is  much  larger  than  the  wave  frequency,  co.  As  a  result 
resonance  only  will  occur  when  the  electrons  and  the  wave  are  traveling  in 
opposite  directions.  Studies  of  many  whistler  spectrograms  show  a  drop  of 
greater  than  10 dB  in  wave  intensity  at  a  frequency  approximately  equal  to  one 
half  the  cyclotron  frequency  of  the  electrons  at  the  magnetic  equator  [16].  VLF 
emissions  are  very  narrow  band  indicating  that  the  magnetic  field  the  electrons 
travel  through  during  resonance  must  be  slowly  varying  over  a  long  length.  For 
these  reasons,  the  interaction  region  is  assumed  to  be  the  section  of  the  field 
aligned  duct  that  is  near  the  magnetic  equator  as  shown  in  Figure  3-5  [17]. 
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Figure  3-5:  Sketch  of  the  paths  taken  by  the  electrons  and  the  wave 
including  the  interaction  region  (not  to  scale). 

If  the  magnetic  field  is  taken  to  be  constant  along  the  interaction  region 
and  £  is  parallel  to  the  z  axis,  equation  3-5  can  be  written  as 

V  =  --2-(Vx£)  ,  3-14 

me 

where  V  is  the  electron  acceleration  expressed  as  the  first  derivative  of  velocity 
with  respect  to  time.  This  vector  equation  can  be  written  in  component  form 
resulting  in  three  differential  equations 

vx  =  -nvy  ,  Vy  =  nvg  ,  v*  =  o  .  3-15 

The  z  component  of  velocity  has  no  time-dependence,  because  the  Lorentz 
force  acts  on  that  part  of  the  electron  velocity  which  is  not  parallel  to  the  field. 
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A  solution  for  these  equations  is 

x  =  rccos(£lt  +  0O)  ,  y  =  rcsin(Qf  +  80)  ,  z  =  -v2t  .  3-16 

The  term  rc  =  vx/a  and  is  the  radius  of  the  cyclotron  orbit,  Vi  is  the  transverse 
electron  velocity,  a  is  the  cyclotron  frequency,  and  e0  is  the  initial  phase  angle. 
Electrons  with  motion  described  by  these  equations  will  be  referred  to  as  the 
"hot"  electrons.  The  resonance  between  these  electrons  and  the  counter 
propagating  wave  will  be  analyzed  through  the  use  of  computer  simulations. 
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IV.  WAVE  PARTICLE  INTERACTION  THEORY 


A.  THE  BACKGROUND  PLASMA 

The  propagation  of  the  VLF  wave  through  the  interaction  region  is  affected 
by  the  properties  of  the  hot  and  cold  plasma  in  the  region.  The  "cold" 
background  plasma  is  more  dense  (pc  =  400/cm3)  than  the  hot  plasma 
described  in  the  previous  chapter.  The  cold  background  plasma  contains  non¬ 
resonance  electrons  that  affect  the  wave  dispersion  relation  while  hot  electrons 
amplify  the  wave.  Without  considering  the  effect  of  the  VLF  wave,  the  cold 
plasma  executes  cyclotron  motion  around  the  magnetic  field  lines.  This  motion 
is  described  by  the  Lorentz  force  equation  and  t  is  given  by 

r=  -^[-cos(flf),  sin(flf),  0]  ,  4-1 

where  r  =  Vci/Q  is  the  radius  of  gyration,  and  the  initial  phase  angle  is  taken  to 
be  zero.  Since  the  cold  electrons  are  non-resonant,  the  small  z  velocity 
component  can  be  ignored.  The  cold  electrons  rotate  in  a  clockwise  direction 
when  looking  along  the  magnetic  field  line  with  a  velocity  equal  to, 

Pc  =  VcJsinfOf),  cos(flf),  0]  .  4-2 


When  the  VLF  wave  is  included,  the  Lorentz  force  equation  is  written  as 


x  {Bw+8)\  , 


4-3 


where  and  0W  are  the  electric  and  magnetic  field  vectors  of  the  propagating 
VLF  wave.  The  magnitude  of  Bw  is  much  less  than  the  geomagnetic  field  in  the 
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interaction  region  and  can  therefore  be  neglected  here.  In  component  form, 
equation  4-3  is 

Vcx  +  nVcy=--%-Ewx  .  4-4 

m 

Vgy  ~  Q ^cjt  =  ~  ^wy  4-5 

Assuming  that  the  incoming  VLF  wave  is  right-hand  circularly  polarized,  the 
fields  are; 

=  Ew[ cos(v),  -  sin(y),  0] ,  Sw  =  S^fsinty),  cos(y),  0]  .  4-6 

where  the  phase  of  the  propagating  wave  is  The  wave  travels 

from  the  south  pole  to  the  north  pole  as  previously  shown  in  Figure  3-5.  The 
wave  propagation  vector  is  taken  to  be  parallel  to  the  geomagnetic  field  and  the 
wave  acts  as  driving  force  for  the  cold  electrons.  As  a  result,  their  steady-state 
motion  must  have  the  same  space  and  time  dependence  as  the  driving  wave. 
If  the  form  for  the  cold  electron  velocity  is  taken  to  be 

Vc  =  Vci[sin(v),  cos(y),  0]  ,  4-7 

and  this  form  is  then  used  in  equations  4-4  and  4-5,  two  equations  which 
include  V'ci  result: 

Vc4to  cos(y)  -  nvci  cos(y)  =  -^Ewcos(y)  ,  4-8 

Vcl(o  sin(y)  -  ClVcl  sin(y)  =  ^E*sin(y)  .  4-9 


Both  equations  can  be  solved  simultaneously  resulting  in  an  expression  for  V.*. 


m{Cb- <b) 


4-10 
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The  current  generated  by  the  motion  of  these  cold  electrons  can  now  be  written 


as 


*  "  m({L)  ts,n't’'  01  • 

4-11 

where  pc  is  the  density  of  the  cold  electrons. 

This  current  density,  along  with  the  expressions  for  Ew 

and  Bw,  can  then 

be  inserted  into  Ampere’s  Law  and  Faraday's  Law 

-  4,1 7  i  1 

VxS„  -  c  -4  *  c  Jf  . 

4-12 

Vxt"~  c  dt  ‘ 

4-13 

Assuming  that  the  VLF  field  envelope,  Ew,  Bw  and  4>,  are  constant  in  time  and 

space,  these  Maxwell  equations  can  be  simplified  to 

o  ^  „  ,  47ie2pc£i* 

B"  E”  m(n-to)  • 

4-14 

Ewck  =  Bwa 

4-15 

The  index  of  refraction  of  the  interaction  region  is  defined  as  n  =  ck/a>.  As  a 

result,  equation  4-15  shows  that  Bw  =  nEw.  Using  this 

substituting  into  equation  4-14  results  in  an  expression  for  n2, 

relationship  and 

©2  ©(ft-©)  ’ 

4-16 

where  the  plasma  frequency  for  the  cold  electrons,  ©c  =  4nezpcim.  The 
expression  for  the  index  of  refraction  given  in  Chapter  II  (equation  2-3)  follows 
directly  from  4-16  because  ©c  »©  and  ©c  Equation  4-16  is  the  exact 
dispersion  relation  for  the  VLF  wave  in  a  cold,  magnetized  plasma  [14]. 
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B.  THE  WAVE  EQUATION 


The  next  step  is  to  develop  a  wave  equation  that  describes  the  evolution 
of  the  VLF  wave  envelope.  It  would  include  the  factors  that  contribute  to  the 
growth  of  the  wave  envelope  as  it  passes  through  the  interaction  region. 
Referring  back  to  the  two  Maxwell  equations  used  in  the  last  section,  the 
current  density  must  now  include  the  hot  electrons 

J=Jc+Jh.  4-17 

where  7h  =  -  eVhl[-  sin(9) ,  cos(6) ,  0]  is  the  current  for  a  single  hot  electron,  Vhi 
is  the  transverse  velocity  of  a  hot  electron  and  0  =  at  +  0O  is  the  phase  of  the 
electron  cyclotron  motion.  The  variation  of  the  wave  envelope  in  space  and 
time  is  assumed  to  be  small  over  one  wavelength  [13], 

Ew  <t  to Ew  ,  <j>  oxt> ,  E'w  «  kEw  ,  «  kty  . 

This  slowly-varying  amplitude  and  phase  approximation  is  appropriate  when  the 
VLF  envelope  describes  a  narrow-band  whistler  mode. 

Returning  to  Ampere’s  Law,  equation  4-12,  and  using  equation  4-6  to 
represent  the  VLF  wave  results  in  the  following  vector  equation: 

8Z£?*[-  cos(\y),sin(\|/),0]+flw(k+3z4>)[sin(\j/),cos(\jr),0]  -  4-18 

— 3fElv[cos(\|f),sin(v),0]+-~(-(i>faf4»)[sin(\|f),cos(v),OJ  =  , 

CO  c 

where  az(..)  =  d(..)/dz  and  3f(..)  =  9(..)/at.  Defining  two  unit  vectors, 

6,  =  [smy  ,cosy  ,0]  ,  e2  =  [cosy  ,-siny  ,0]  ,  4-19 

and  projecting  Ampere’s  Law  onto  these  orthogonal  unit  vectors  results  in  two 
scalar  first-order  differentia’  equations, 
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(D^fr 

ckBw  -  (oEw  +  cBwdz<b  +  -  47ie1/’/,4.cos(\|f+8)  ,  4-20 

•fcu) 

cdzBw  +  dtEw  =  -4jieVh.|Sin(\|/+8)  .  4-21 

Using  nEw  -  Bw  and  the  dispersion  relation,  4-16,  further  simplifies  these 
equations  to  obtain, 

ncBwdzt (>  +  Bwd,$  =  -4nenVhlcos(y+Q)  ,  4-22 

ncdzBw  +  d,Bw  =  -4nensin(y+B)  .  4-23 

At  this  point  it  is  convenient  to  define  the  quantity  b : 

b  =  -^r-exp(/4>)  .  4-24 

me 

The  complex  field  b  has  units  of  inverse  seconds  (s-1),  and  is  =  1/s,  when  the 
VLF  wave  amplitude  is  at  a  maximum  (lOpT).  It  is  anticipated  that  the  VLF 
fields  will  not  be  bigger  than  this  value.  Taking  the  partial  derivative  of  b  with 
respect  to  z  and  t  in  accordance  with  the  chain  rule  results  in 

dzb  =  —exp{i$)[dzBw  +  iBwdz§]  ,  4-25 

me 

Btb  =  —  exp(/»[3fB*  +  iBwdz$]  .  4-26 

me 

These  expressions  can  then  be  used  to  organize  4-22  and  4-23  into  a  complex 
equation, 

(ncdz  +  dt)b  =  Vexp[-/(C+e)] ,  4-27 

where  =  kz  -  at  +  at  is  the  longitudinal  electron  phase  and  V=nVhl/c  is  the 
dimensionless  transverse  electron  velocity. 
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The  current  in  Ampere’s  Law  is  for  a  single  hot  electron  in  the  presence  of 
the  incoming  VLF  wave.  The  current  for  a  beam  of  hot  electrons  is  obtained  by 
averaging  over  the  transverse  velocity  V,  the  transverse  phase  e,  and  the 
longitudinal  phase  £•  The  average  is  symbolized  by  <...>,  and  must  be  weighted 
by  the  hot  electron  beam  density  p„.  The  wave  equation  4-27  for  a  beam  of 
hot  electrons  is  then  as 

( ncdz  +  d,)b  =  -/ co2/, <  Vexp[-/ (C+0)]>  .  4-28 

where  a2h  =  4ne2ph/m  is  the  hot  electron  plasma  frequency.  When  the  VLF 
wave  envelope  is  spatially  flat  over  many  wavelengths,  d2b  is  zero  and  the 
wave  equation  can  be  written  as 

b  =  -/©V  Vexp[-/(C+0)>  .  4-29 

As  the  incoming  VLF  wave  enters  the  interaction  region  it  encounters  a 
random  distribution  of  hot  electrons  whose  z  motion  is  described  by 
zh  =  Zho  ~  vhzt ■  The  current  contains  the  factor  <exp[-i(kz0)]>  where  z0  is  the 
initial  random  position  of  a  hot  electron.  The  number  of  hot  electrons  in  the 
average  <...>  describing  one  wavelength  of  the  incoming  wave  is  =  1013  so  that 
the  current  will  average  to  a  value  very  close  to  zero.  The  right  side  of 
equation  4-29  is  small  regardless  of  the  value  of  co/,  and  therefore,  the  VLF 
wave  is  not  amplified.  What  is  needed  is  a  description  of  how  the  positions  of 
the  hot  electrons  evolve  as  the  wave  passes.  The  dynamical  equations  for  the 
hot  electrons  in  the  presence  of  the  VLF  wave  are  derived  next. 
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C.  DEVELOPMENT  OF  ELECTRON  EQUATIONS  OF  MOTION 

The  motion  of  the  hot  electrons  in  the  presence  of  the  incoming  field,  like 
the  motion  of  the  cold  electrons,  is  governed  by  the  Lorentz  force  equation 


=  -  *—•  [£*  +  Vhx(&w  +  £)] 
m 


4-30 


but  unlike  the  analysis  used  for  the  cold  electrons,  the  magnetic  field  of  the 
incoming  wave  Bw  is  not  dropped.  The  geomagnetic  field  is  assumed  to  be  of 
the  form  £?  =  [0,0,6],  and  Ew  and  Bw  are  taken  from  equation  4-6,  so  that 
equation  4-30  becomes 


cosy  + 


4-31 


r  i 

Vhy  =  — 
y  m 

-Ew 

c|  o 
_ * 

siny-  ^Vhx 

4-32 


Ew^(Vhx cosy  -  V'tySiny) 


4-33 


where  y  =  /cz-mf-Hji.  The  hot  electron  velocity  in  the  x-y  plane  can  be  expressed 
in  polar  coordinates  in  terms  of  and  eA, 

Vhx=-Vhisin(ei)  ,  Vhy  =  V^cosiBi)  .  4-34 

Expressions  for  Vhx  and  Vhy  can  be  obtained  from  equation  4-34  through  the 
use  of  the  chain  rule.  After  making  the  appropriate  substitutions  into  equations 
4-31  and  4-32  the  following  expressions  result: 

Vi hz  =  ^~nVhl  sin^  +  y)  ,  4-35 

Vtn  =  ~~(c-Vhzn)s\n(Ql  +  y)  ,  4-36 
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=  Q  +  — if-  val  {c-V^n  )cos(ei  +  v)  .  4-37 

me 

Making  the  definition  0^  =  cit  +  6  the  argument  of  the  sin  and  cos  functions 
in  equations  4-35,  4-36  and  4-37  becomes 

0^  +  v  =  kzh  -  (at  +  Clt  +  0  +  <|>  .  4-38 

The  electron  phase  is  C  =  te/,  -  cof  +  Clt.  Using  0A  =  n  +  0,  the  hot  electron 
equations  of  motion  take  the  following  form; 

'ih  =  VAAsin(C  +  0  +  <|»)  ,  4-39 

me 

r  ’* 

1--^sin(C  +  e  +  «  .  4-40 

nm  c 

V. 

r  *■ 

8  =  7^77“  1  "  ~~T~  cos(C  +  8  +  4>)  .  4-41 

nmVhl  c  J 

where  Bw  =  nEw  and  zh  =  Vh2. 

Looking  at  the  definition  for  it  is  clear  that  £  =  #cz.  Making  this 
substitution  into  equation  4-39  results  in 

C  =  ^-co-^sintC  +  0  +  0)  .  4-42 

me  c 

Equation  4-42  can  be  further  simplified  using  the  definitions  for  a  complex  field 
b  =  Ible'*  where  lb  I  =  eBjmc  and  where  V  =  nVhi/c.  The  expression  for  the 
electron  phase  is  now 

£=  Va>lblsin(C  +  0  +  0)  .  4-43 

Equation  4-43  will  be  referred  to  as  the  pendulum  equation  for  the  electrons. 

Equations  4-40  and  4-41  can  also  be  simplified  through  the  use  of  the 
phase  and  group  velocity  of  the  VLF  wave 
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4-44 


i/  <i  (a  y  _  dc o 
p  =  n  =  k  '  dk 

where  ck  =  n to  and  n  is  the  index  of  refraction  in  the  interaction  region.  Using 
the  index  for  the  case  of  o)c:xo,n, 


-  oxflc 
Vco(fi-<o) 

the  group  velocity  can  be  evaluated  in  terms  of  n  as 


4-45 


4-46 


The  phenomenon  of  a  "cutoff"  frequency  was  described  in  Chapter  III,  and 
is  defined  as  one  half  the  electron  cyclotron  frequency,  n/2,  at  the  magnetic 
equator.  The  simulations  presented  in  the  next  chapter  will  assume  that  the 
frequency  of  the  incoming  VLF  wave  is  to  =  n/2.  When  this  is  substituted  into 
equation  4-46,  it  is  found  that  Vg  ~  Vp.  At  resonance  and  at  cutoff,  =  -Vp  so 
that  the  electrons  and  the  wave  travel  at  equal  speeds  in  opposite  directions. 
Using  -cm  for  Vh2  in  equations  4-40  and  4-41  as  well  as  the  definitions  for  the 
magnitude  of  lb  I  and  V  results  in 


V  =  2lblsin(£  +  0  +  <|>)  , 

4-47 

e  =  -^-cos(C  +  e  +  <t>)  . 

4-48 

The  equations  of  motion  for  the  electron  in  the  presence  of  the  VLF  wave 
equations,  4-43,  4-47,  4-48,  along  with  the  wave  equation,  4-29,  are  now  in  a 
convenient  form  to  integrate  numerically  and  understand  further. 

A  final  modification  to  the  equation  for '(  comes  when  the  electron  phase 
velocity  is  defined  as  v  =  £.  This  allows  all  the  important  equations  describing 
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the  hot  electrons  and  the  VLF  field  to  be  in  the  form  of  first-order  differential 
equations.  The  collected  set  of  equations  are 


v  =  V<a\b  lsin(£  +  0  +  $)  , 

4-49 

C  =  v  . 

4-50 

V  =  2lblsin(£  +  0  +  4>)  , 

4-51 

e  =  ^-cos(C  +  e  +  <i>)  , 

4-52 

6  =  -/  co/,  2<  Ve~'^  +  e)>  . 

4-53 
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V.  SIMULATIONS  AND  DISCUSSION 


A.  EXPLANATION  OF  NUMERICAL  INTEGRATION  METHODS 

The  equations  presented  at  the  end  of  Chapter  IV  are  non-linear,  so  that 
the  only  method  of  solution  is  through  the  use  of  numerical  techniques.  This 
process  can  be  carried  out  on  a  desk-  top  computer.  Assuming  a  small  finite 
time-step.  At,  the  electron  equations  can  be  written  as 


Av  =  Vcol  b  lsin(£  +  0  +  $)A  t  , 

5-1 

AC  =  vA t  , 

5-2 

A  V  =  210  lsin(C  +  0  +  4>)Af  , 

5-3 

A0  =  ^-cos(C  +  0  +  , 

5-4 

The  equations  for  Av  and  a;  each  have  a  factor  of  u  which  is  not  present  in  the 
equations  for  AV  or  A0.  The  value  of  co  is  ~104/s,  and  it  is  therefore 
anticipated  that  the  equation  for  v  will  evolve  much  faster  than  the  last  two 
equations.  As  a  result,  the  first  two  equations  will  be  numerically  integrated 
using  the  fourth-order  Runge-Kutta  method  [18]  while  the  simpler,  more  efficient 
Euler-Cromer  method  [18]  will  be  used  to  integrate  the  last  two  equations. 
Each  equation  describes  how  the  values  for  each  variable  change  over  the 
short  time  interval  At.  The  simulation  will  use  this  process  to  keep  track  of  and 
plot  how  the  electrons  evolve  over  time. 

The  wave  equation  is  also  a  first-order  differential  equation  in  time,  but  V,  £ 
and  6  on  the  right-hand  side  of  4-53  must  be  averaged  over  all  the  sampled  hot 
electrons.  In  the  simulation,  the  values  for  v,  C,  V  and  9  will  be  found  first  for 
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each  of  the  hot  electrons.  The  average  values  for  V,  ;  and  e  will  then  be 
determined,  and  used  in  finding  a  new  value  for  b.  For  a  small  finite  time-step 
At,  the  wave  equation  can  be  written  in  complex  form  as 

A b  =  -i(oh  2<  I'e-'K  +  ®>>A t  .  5-5 

The  computer  actually  evaluates  this  equation  as  written  in  real  and  imaginary 
components.  The  actual  equations  used  by  the  computer  are 

A br  =  -tO/,2<Vsin(£  +  0  +  4>)>Af  ,  5-6 

Abj  =  -<d,,2<Vcos(£  +  0  +  $)>At  .  5-7 

where  b  =  br  +  ibr  Using  these  real  and  imaginary  components  in  the 


equations  for  the  electron  motion  results  in 

Av  =  Vco[b,sin(C  +  0)  +  b,cos(£  +  0)]Af  ,  5-8 

A^  =  vA  t  ,  5-9 

A  V  =  2[brsin(C  +  0)  +  b,  cos(C  +  0)]Af  ,  5-10 

A0  =  -^[b,cos((;  +  0)  -  b(sin(C  +  0)]Af  .  5-1 1 


These  are  the  equations  which  will  be  integrated  by  the  computer  using  the 
methods  mentioned  above. 

B.  INITIAL  CONDITIONS 

The  position  of  each  electron  and  the  amplitude  and  phase  of  the  VLF 
wave  can  now  be  determined,  but  the  initial  values  for  each  of  the  variables  is 
required  for  integration.  As  stated  in  Chapter  IV,  there  are  many  electrons 
(=  1013)  in  a  volume  of  the  interaction  region  one  wave  length  long.  It  is  not 
possible  to  integrate  this  many  electrons  so  that  a  smaller  sampling  will  be 
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used.  These  electrons  are  assumed  to  be  in  a  beam  with  a  uniform  distribution 
in  initial  C  and  initial  6  ranging  between  -n  to  ic.  The  initial  values  for  v  will  be 
determined  by  the  distribution  functions  f  (v),  where  ov  characterizes  the  spread 
in  f  (v),  while  v0  characterizes  the  function  f(v).  The  distribution  function  g{V) 
provides  the  initial  distribution  in  the  dimensionless  transverse  velocity  V  with 
av  characterizing  the  spread  in  g(V),  and  V0  characterizes  the  position  of  g(V). 
Since  the  equation  of  motion  for  e  involves  V~\  it  is  important  that  g(V) 
approach  zero  for  small  V.  The  loss  cone  for  mirroring  electrons  in  the 
ionosphere  assures  this  property. 

The  initial  VLF  field  is  b,  =  b0  and  b,  =  0  so  that  the  initial  phase  angle  is 
<|>(0)  =  0.  The  plasma  frequency,  <ah,  depends  on  the  density  of  hot  electrons  at 
the  time  of  the  VLF  interaction  in  the  Earth's  ionosphere  and  can  vary 
significantly  with  the  level  of  magnetic  activity.  When  ph  =  i/m3,  <oh  =  50/s  [19] 
provides  a  typical  value  for  the  ionosphere,  but  other  values  will  be  examined  in 
the  simulations.  The  length  of  time,  T,  over  which  the  interaction  takes  place  is 
typically  a  fraction  of  a  second,  based  on  the  observation  that  whistler  traces 
are  approximately  one  second  long.  The  interaction  occurs  in  a  relatively  small 
portion  of  the  total  length  of  a  field  aligned  duct,  therefore  the  run  time  for  a 
simulation  should  be  a  fraction  of  that  total  time. 

C.  INTERACTIONS  FOR  MONO-ENERGETIC  BEAM 

The  first  set  of  simulations  are  for  the  case  were  f  (v)  and  g{V)  are  delta- 
functions,  and  all  the  electrons  start  with  the  same  initial  value  for  v0  =  0  and 
VQ  =  1  corresponding  to  a  pitch  angle  of  n/4.  Figure  5-1  shows  the  results  of 
this  simulation  for  a  VLF  wave  with  a  frequency  ©=  104/s.  The  plot  on  the  left 
shows  the  electron  phase  velocity,  v,  versus  the  electron  phase,  £.  In  the 
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Figure  5-1 :  Electron  phase-space,  gain  and  phase  for  weak  fields  and 
low  gain. 


phase-space  plot,  the  path  and  final  positions  for  each  of  the  electrons  are 
plotted.  Those  paths  which  move  down  in  v  indicate  that  these  electrons  are 
giving  up  energy  to  the  wave,  while  those  moving  up  in  v  are  increasing  in 
energy,  removing  it  from  the  wave.  In  this  case,  almost  as  many  electrons  go 
up  in  energy  as  go  down.  The  electrons  clearly  bunch  around  two  values  of  ; 
in  one  wavelenght  of  the  VLF  wave.  The  right-hand  plots  show  how  the  gain, 
6,  and  phase  of  the  VLF  wave  ,  $,  develop  over  the  course  of  the  interaction. 
The  gain  is  a  measure  of  how  the  wave  amplitude  changes  over  time  and  is 
defined  as 


5-12 


The  average  on  the  right  side  of  the  wave  equation,  4-53,  provides  the  driving 
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current  for  the  VLF  wave.  When  there  is  bunching  the  average  is  close  to  zero, 
and  the  wave  does  not  grow  significantly.  When  bunching  finally  occurs  near 
the  end  of  the  undulator,  the  gain  increases  to  a  final  value  of  G  =  3%.  The  size 
of  lb  I  is  small  and  does  not  increase  significantly  because  <oh  is  small.  The 
phase  of  the  VLF  wave,  4>,  does  not  change  significantly  during  the  interaction. 

In  Figure  5-2,  the  same  initial  conditions  are  used  but  the  plots  of  G  and  4 
are  replaced  by  a  plot  showing  the  evolution  of  V  and  e  on  the  right. 
Throughout  the  integration,  these  variables  remain  very  close  to  their  initial 
values. 


Figure  5-2:  Electron  phase-spaces  (O)  and  (6,V)  for  weak  fields  and 
low  gain. 

In  Figure  5-3  the  hot  plasma  frequency  is  increased  to  a*  =  5/s  resulting  in 
an  increase  in  gain  for  this  interaction.  At  the  start  of  the  interaction  the 
electrons  begin  to  bunch  as  in  the  previous  cases.  This  causes  +  to  increase 
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so  that  the  point  around  which  the  electrons  bunch  advances  a  little  along  £• 
The  larger  driving  current  increases  the  amplitude  of  the  VLF  wave  and  the 
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Figure  5-3:  Electron  phase-space,  gain  and  phase  for  weak  fields  and 
high  gain. 


gain.  The  gain  increases  exponentially  in  time  and  reaches  a  final  value  of 
28 dB.  The  phase  increases  significantly  to  nearly  n  at  the  end  of  the 
interaction.  The  final  size  of  lb  I  is  still  small  but  it  is  much  greater  than  the  initial 
value  due  to  high  gain.  This  set  of  conditions  produced  high  gain  in  relatively 
weak  fields. 

Figure  5-4  shows  that  again  the  variables  V  and  6  did  not  evolve 
significantly. 
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Figure  5-4:  Electron  phase-spaces  (C,v)  and  (Q.V)  for  weak  fields  and 
high  gain. 

The  next  set  of  simulations  study  the  effect  of  large  VLF  field  strengths.  In 
Figure  5-5,  the  value  of  b0  =  0.1  /s  is  increased  by  a  factor  of  ten  over  the  value 
in  Figure  5-1  while  o*h  =  0.5/s  stays  the  same.  This  larger  initial  field  causes  the 
electrons  to  over-bunch  as  shown  in  the  electron  phase-space  plot  on  the  left. 
As  the  bunching  starts,  the  gain  increases,  but  as  the  over-bunching  continues 
the  rate  of  increase  in  the  gain  slows  and  finally  levels  off  at  the  end  of  the 
simulation.  As  the  electrons  move  past  their  bunched  positions,  they  take 
energy  away  from  the  VLF  wave  reducing  the  amplitude.  The  final  gain  is  only 
G  =  0.4%  for  the  same  value  of  a>h  as  in  Figure  5-1 ,  but  it  has  decreased 
because  of  large  Ibl.  This  also  leads  to  almost  no  change  in 
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Figure  5-5:  Electron  phase-space,  gain  and  phase  for  strong  fields 
and  low  gain. 

In  the  phase-space  plot  of  Figure  5-5,  the  boundary  lines  define  the 
maximum  radius  within  the  closed  orbit  region.  The  boundary  lines  are  located 
at  Av  =  ±2Vlbi(oV  above  and  below  resonance.  The  frequency  at  which  the 
electrons  orbit  about  fixed  points  within  the  closed  orbit  region  is  known  as  the 
"synchrotron  frequency",  vs  =  VIbItoV.  It  is  possible  that  the  gain  can  also 
oscillate  at  this  synchrotron  frequency. 

In  Figure  5-6  the  plot  of  electron  velocity  versus  6  does  not  evolve  during 
the  interaction. 
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Figure  5-6:  Electron  phase-spaces  (C.v)  and  (Q,V)  for  strong  fields  and 
low  gain. 

Figure  5-7  shows  the  results  of  high  gain  with  (awubh  =  5/s  and  a  strong 
field  b0  =  0.3/s.  The  electrons  over-bunch  again,  but  the  gain  is  much  higher 
because  a>h  is  bigger.  The  over-bunching  limits  the  gain  to  G  =  5.9dB 
compared  to  G  =  29 dB  in  Figures  5-3  and  5-4.  Some  of  the  electrons  crossed 
the  lower  boundary  of  the  closed-orbit  region  because  the  value  of  lb  I  is  larger 
at  f=0.06s  than  the  final  value  used  to  draw  the  boundary  lines  at  ±2^\b\aV. 
The  gain  plot  shows  a  tendency  for  the  gain  to  oscillate  at  the  synchrotron 
frequency.  The  phase  change,  is  large  as  expected  with  a  large  value  of  a>„. 
The  data  collected  on  VLF  waves  in  the  magnetosphere  indicates  the 
interaction  involves  strong  fields  with  high-gain.  In  Figure  5-8,  the  plot  of 
electron  transverse  velocity  V  versus  e  still  does  not  evolve  during  the 
interaction. 
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Figure  5-7:  Electron  phase-space,  gain  and  phase  for  strong  fields 
and  high  gain. 
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Figure  5-8:  Electron  phase-spaces  (C,v)  and  (e,V)  for  strong  fields  and 
high  gain. 

Figure  5-9  shows  the  same  strong-field,  high-gain  simulation  as  in  Figure 
5-7,  but  the  electrons  are  moved  off-resonance  to  v0=  130/s.  The  electrons 


follow  open  orbit  paths  which  are  outside  the  closed  orbits  region.  The  large 
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initial  value  for  v0  results  in  a  large  initial  £  where  C  =  Co  +  vor-  This  large  £ 
oscillates  so  quickly  that  the  average  of  £  is  zero  over  one  phase  $.  This 
results  in  weak  coupling  to  the  VLF  field  and  therefore  the  gain  in  this  case  is 
smaller,  G  =  3.800,  than  in  Figure  5-7  where  G  =  5.900.  The  plot  of  V  versus  9 
in  Figure  5-10  does  not  evolve  during  the  interaction. 
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Figure  5-9:  Electron  phase-space,  gain  and  phase  for  strong  fields 
and  high  gain  starting  off-resonance. 
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Figure  5-10:  Electron  phase-spaces  (C,v)  and  (e,V)  for  strong  fields 
and  high  gain  starting  off-resonance. 


D.  SIMULATIONS  USING  DISTRIBUTIONS 

The  electrons  in  the  magnetosphere  are  not  monoenergetic.  Studies  of 
data  obtained  from  spacecraft  flying  through  the  magnetosphere  indicate  that 
the  density  of  electrons  between  Eh  =  750eV  and  Eh  =  50keV  varies  as  Eh~n 
where  n  can  have  values  from  1.5  to  3.5  [20].  There  are  many  more  electrons 
with  low  energies  than  there  are  at  higher  energies  for  this  type  of  distribution 
function.  Using  the  relationship  between  energy  and  velocity  means  the 
number  density  of  electrons  in  terms  of  velocity  is  Vh~2n.  This  function  can  also 
represent  the  distribution  of  hot  electrons  in  terms  of  and  because  these 
quantities  are  directly  related  to  the  electron  velocity.  The  initial  phase  velocity 
of  the  electrons,  v0,  is  directly  related  to  the  initial  velocity  in  the  z  direction,  Vhz. 
When  the  electrons  have  a  value  of  v0  which  is  significantly  away  from 
resonance,  there  is  poor  coupling  to  the  VLF  field,  and  the  gain  is  small  as 
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shown  in  Figure  5-9.  As  a  result  only  those  electrons  with  v0  in  a  narrow  band 
around  resonance  need  be  considered.  The  size  of  this  narrow  band  around  v0 
is  related  to  the  spread  in  Vh2  through  the  relationship, 


Avq  = 
to  cm 


5-16 


where  cm  is  the  velocity  of  resonant  electrons.  When  Av  =  500/s  and 
(o  =  15000/s  the  spread  around  the  resonant  velocity  is  only  ±3%.  This  is  such  a 
small  portion  of  the  overall  distribution,  Eh~n,  that  the  initial  electron  phase 
velocity  distribution,  /(v0),  can  be  represented  by  a  linear  function  with  a 
negative  slope.  Smaller  values  of  v0  will  have  a  higher  probability,  because 
lower  velocities  are  more  likely  than  higher  velocities  for  this  type  of  distribution 
function.  The  actual  f(v0)  used  is  triangular  in  shape  with  a  full-width-at  half¬ 
maximum,  ov,  wider  than  the  window  of  v0  values  used.  The  result  is  a  sloped 
distribution  in  /(v0)  used  to  start  the  simulation  electrons. 

The  distribution  function  g(V)  used  to  select  the  initial  values  for  V  must 
reflect  the  properties  of  the  velocity  distribution  discussed  above  as  well  as  the 
distribution  in  pitch  angles  a  because  V  =  nVhl/c.  The  lower  boundary  for  g{V) 
is  defined  by  the  loss  cone,  am.  Electrons  with  pitch  angles  less  than  am 
precipitate  out  of  the  field  aligned  ducts  and  are  not  present  to  contribute  to 
wave  growth.  Because  this  angle  is  so  small,  Vhl~  Vh  and  Vhz  =  Vh.  The 
resonant  velocity  is  defined  as  cm,  therefore,  at  am 


V  = 


cm 


a, 


m 


5-17 


An  exponential  distribution  will  be  used  for  g(V)  with  the  peak  at  V0  which 
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corresponds  to  a  =  am.  The  full  width  at  half  maximum,  av  =  1,  corresponds  to 
a  pitch  angle  a  =  45°,  where  Vhi  = 

In  Figure  5-11,  the  simulation  is  run  for  initial  parameters  which  are 
intended  to  represent  the  environment  in  the  interaction  region.  The  initial 
values  for  v  and  V  for  each  electron  is  determined  using  f(v)  and  g{V).  The 
value  the  hot  electron  plasma  frequency  is  taken  to  be  <ah  =  25/s  corresponding 
to  a  density  of  hot  electrons,  p/,  =  .18/m3  in  the  magnetosphere.  The  initial  VLF 
field  is  b0  =  .01  /s  which  is  lower  than  the  strong  field  cases  previously  presented 
in  Figures  5-7  through  5-10.  The  spread  in  Vz  is  found  using  5-16  and  is  equal 
to  ±3%  of  the  resonant  velocity.  A  large  number  of  particles  (8000)  are  used  in 
the  simulations,  therefore  only  the  final  positions  are  plotted  in  phase-space. 
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Figure  5-11:  Electron  phase-space  gain  and  phase  for 
magnetospheric  initial  conditions,  short  interaction  time. 
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In  phase-space,  particles  with  lower  energy  are  at  the  bottom.  In  accordance 
with  the  description  of  / (v)  and  g(V),  there  are  more  electrons  in  the  lower  half 
or  the  phase-space  plot  than  in  the  upper  half.  Although  this  may  be  difficult  to 
perceive  from  just  looking  at  the  plot,  the  simulation  checks  for  this  by  counting 
the  number  of  electrons  in  a  small  strip  at  the  top  and  bottom  of  phase-space 
and  taking  the  ratio  of  these  numbers.  In  Figure  5-11,  the  gain  grows 
exponentially  and  the  phase  does  not  change  significantly. 

When  the  simulation  runs  for  a  longer  integration  time,  T,  as  in  Figure  5- 
1 2,  the  gain  saturates  then  begins  to  oscillate.  This  shows  the  saturation  of  the 
VLF  wave  in  strong  VLF  fields. 
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Figure  5-12:  Electron  phase-space  gain  and  phase  for  magnetospheric  initial 
conditions,  long  interaction  time. 
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VI.  SUGGESTIONS  FOR  FUTURE  WORK 


A.  IMPROVEMENTS  TO  THE  SIMULATION 


The  simulation  at  this  point  in  development  assumes  the  geomagnetic 
field  is  homogeneous  throughout  the  interaction  region.  A  better  approximation 
would  include  a  factor  in  the  electron  equations  of  motion  which  accounts  for 
the  increasing  field  seen  by  the  electrons  as  they  move  to  higher  latitudes.  This 
increase  in  the  field  decelerates  the  electrons  in  the  z  direction  and  changes 
the  VLF  resonance  condition.  The  expression  for  this  deceleration  was 
presented  in  Chapter  III, 


Vhi2  dB(z) 
28  dz 


6-1 


where  8(z)  is  the  geomagnetic  field  as  a  function  of  position  z  [15J.  In  the 
absence  of  the  VLF  interaction,  the  velocity  of  the  hot  electrons  at  the  magnetic 
equator  is  constant  in  time,  and  can  be  expressed  in  component  form  as 

6-2 


Taking  the  time  derivative  of  both  sides  in  6-2,  and  using  the  expression  for 
z  =  <  results  in  an  expression  for  the  transverse  components 

vh,  =  aBW  .  6-3 

*J  20  dz 

These  factors  can  now  be  added  to  the  expressions  found  previously  for  zh  and 
V'/.i  giving, 
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6-4 


Zh 


eBw 

me 


V^sinfC  +  e  +  4>)  - 


V*i  96  (z) 
2  8  dz 


eBw 

nm 


sin(£  +  6  +  $)  + 


VhzVhj  dB(z) 
26  dz 
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The  field  around  the  equator  as  a  function  of  z  is 


S(z)  =  B0 


1  + 


9  z2 


12  flL 
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where  Re  is  the  radius  of  the  earth  [21],  The  development  from  here  on  is  the 
same  as  previously  presented  in  Chapter  IV. 

This  addition  will  cause  the  cyclotron  frequency  of  the  hot  electrons  to 
change  as  the  electrons  travel  along  their  helical  paths.  As  a  result,  the 
electrons  move  in  and  out  of  resonance  with  the  wave.  When  the  field  is 
assumed  to  be  homogeneous  the  interaction  turns  on  and  off  abruptly.  By 
including  the  effect  of  the  changing  field  the  interaction  turns  on  and  off 
gradually  in  a  similar  way  to  that  in  the  magnetosphere. 

The  distribution  functions  /(v)  and  g{V)  present  another  area  which  needs 
further  investigation.  An  accurate  description  of  the  initial  energies  of  the 
resonant  electrons  will  improve  the  output  of  the  simulation.  This  problem  is 
complicated  because  these  functions  vary  significantly  with  magnetic  activity.  If 
all  other  initial  conditions  were  accurately  known,  various  types  of  distribution 
functions  could  be  studied  and  then  compared  to  experimental  data  in  hopes  of 
predicting  the  electron  density. 

In  all  the  cases  examined  in  Chapter  V,  the  equations  for  V  and  6  did  not 
evolve.  Future  simulations  may  not  need  to  integrate  these  equations  and  still 
produce  meaningful  results.  In  that  case,  this  type  of  simulation  would  closely 


53 


resemble  those  used  in  the  study  of  Free  Electron  Lasers  (FEL).  Each  would 
integrate  a  pendulum  equation,  following  the  development  of  the  electron  phase 
velocity  and  the  coupling  to  the  field.  Even  if  the  simulation  integrates  the 
single  equation,  the  distribution  function  g(V)  must  still  be  included.  In  the  FEL, 
all  electrons  have  the  same  Vi  which  is  directly  related  to  V,  but  that  is  not  the 
case  in  the  case  of  VLF  interactions  in  the  magnetosphere. 

B.  COMPARISON  TO  FREE  ELECTRON  LASERS 

The  theory  and  simulations  presented  up  to  this  point  are  all  in  an  attempt 
to  describe  what  is  a  naturally  occurring  effect.  The  process  used  to  amplify 
light  in  the  Free  Electron  Laser  (FEL)  is  very  similar  to  the  amplifying  process 
described  for  VLF  waves  in  the  magnetosphere.  Any  parallels  in  theory  are 
beneficial  because  an  FEL  is  a  laboratory  device  making  it  much  easier  to  study 
topics  of  mutual  interest.  The  areas  common  to  both  subjects  include: 

1 .  The  FEL  amplifier  couples  the  radiation  coming  from  electrons  accelerated 
in  helical  orbits  to  a  transverse  laser  light  field.  The  electrons  are 
accelerated  as  a  result  of  a  transverse  magnetic  field,  but  the  kind  of 
paths  they  follow  is  the  same  as  for  electrons  in  the  magnetosphere. 

2.  The  radiation  produced  by  the  FEL  is  very  narrow  band  and  can  result  in 
a  gain  of  several  thousand. 

3.  The  evolution  of  the  electron  phase-space  can  be  described  using  the 
pendulum  equation:  v  =  la  lcos(C  +  $),  where  (  =  v.  Simulations  used  to 
investigate  aspects  of  FEL  theory  output  the  same  type  of  displays  as 
used  in  Chapter  V. 
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4.  An  FEL  operating  in  the  microwave  region  makes  use  of  a  waveguide 
resulting  in  a  dispersion  relation,  run  =  kc.[ 22] 

5.  The  electron  distribution  function  is  important  to  FEL  and  the  VLF 
interaction. 

6.  When  the  transverse  field  in  the  FEL  is  tapered  at  the  ends  of  the 
interaction  length,  the  effect  is  similar  to  that  when  the  magnetic  field 
varies  in  the  magnetosphere. 

7.  Advanced  effects  such  as  saturation,  particle  trapping  and  non-linear 
dynamics  are  all  common  to  both  systems. 
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VII.  CONCLUSIONS 


An  overall  description  of  whistler  waves  and  the  related  topic  of  VLF 
emissions  is  presented  in  Chapter  II  along  with  a  discussion  of  how  these 
phenomena  can  change  the  properties  of  the  ionosphere  and  thereby  degrade 
the  performance  of  communication  and  navigation  systems  used  by  the  U.  S. 
Navy.  A  method  using  whistlers  to  remotely  determine  the  electron  density  in 
the  plasmasphere  is  also  presented  here.  The  theory  of  wave-particle 
interaction  is  presented  to  explain  how  the  waves  are  amplified.  The 
propagation  properties  of  the  magnetosphere  are  presented  resulting  in  a 
dispersion  relation  for  waves  traveling  through  a  cold,  magnetized  background 
plasma.  The  equations  of  motion  for  electrons  making  helical  paths  around 
magnetic  field  lines  are  presented  at  the  end  of  Chapter  III. 

A  self-consistent  wave  equation  for  VLF  signals  passing  through  a  cold, 
magnetized  plasma  is  developed  in  Chapter  IV.  A  new  set  of  equations  for  the 
resonant  electrons  in  the  presence  of  the  VLF  field  is  also  presented  at  the  end 
of  the  chapter.  This  set  of  nonlinear  differential  equations  are  solved  using 
numerical  integration  carried  out  on  a  desk  top  computer.  The  first  set  of 
simulations  are  for  a  monoenergetic  beam  of  resonant  electrons.  Weak  and 
strong  initial  field  amplitudes  are  examined  for  both  high  and  low  gain  cases. 
The  outputs  show  the  electron  phase-space  paths  as  well  as  plot  of  gain  and 
the  phase  of  the  field  as  a  function  of  time.  For  weak  fields,  the  electrons 
bunch  around  two  points  in  one  wave  length  of  the  VLF  signal.  The  amount  of 
gain  is  governed  by  the  density  of  the  hot  electrons.  In  strong  fields,  the 
electrons  over-bunch,  taking  energy  away  from  the  VLF  field.  The  gain 
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saturates,  then  oscillates  at  the  synchrotron  frequency.  Distribution  functions  to 
represent  the  density  of  the  hot  electrons  in  the  magnetosphere  are  used  to 
provide  a  more  realistic  set  of  simulations.  A  large  number  of  electrons  are 
needed  so  only  the  final  positions  are  plotted  in  phase-space.  The  gain  grows 
exponentially,  reaches  saturation  and  then  begins  to  oscillate. 

A  number  of  suggestions  are  presented  in  Chapter  VI  to  make  the 
simulation  more  realistic.  Many  aspects  of  Free  Electron  Laser  theory  are  very 
similar  to  that  presented  for  wave-particle  interactions  in  the  magnetosphere. 
These  similarities  are  listed  at  the  end  of  Chapter  VI. 
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